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Abstract 

An analytical model previously developed to study the structure of the 
magnetic field for the TEXTOR-DED [S.S. Abdullaev et al. Phys. Plas- 
mas, 6, 153 (1999)] is applied to the similar study of the Ergodic Divertor 
of Tore Supra tokamak [Ph. Ghendrih, Plasma Phys. Control. Fusion, 38, 
1653 (1996)]. The coil configuration of ED Tore Supra consists of six mod- 
ules equidistantly located along the toroidal direction on the low-field-side 
of the torus with given toroidal and poloidal extensions. The Hamiltonian 
formulation of field line equations in straight-field-line coordinates (Boozer 
coordinates) and the computationally efficient mapping method for integra- 
tion of the Hamiltonian field line equations are used to study the magnetic 
field structure in the ED. Asymptotical formulas for the perturbation mag- 
netic field created by the ED coils are obtained and the spectrum of magnetic 
perturbations is analyzed and compared with the one of the TEXTOR-DED. 
The structure of ergodic and laminar zones are studied by plotting Poincare 
sections, so-called laminar plots (contour plots of wall to wall connection 
lengths) and magnetic footprints. The radial profiles of field line diffusion co- 
efficients are calculated for different perturbation currents and it is found that 
for the Tore Supra case in the ergodic zone the numerical field line diffusion 
coefficients perfectly follow the quasilinear formula for smaller perturbation 
currents although the situation is different for the maximum perturbation cur- 
rent. 
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1 Introduction 



The problem of plasma-wall interaction is one of the key issues in nuclear fusion 
research. To control the plasma edge, the concept of an ergodic divertor has been 
proposed in Refs. [1-3]). The idea of the ergodic divertor is based on the creation 
of a stochastic zone of magnetic field lines open to the plasma wall which would 
guide the ionized particles and energy to special divertor plates, and prevent the 
penetration of wall-released impurities into the plasma core. The stochastic field 
lines at the plasma edge should be formed by externally created resonant magnetic 
perturbations using external coils. 

The ergodic divertor (ED) (ergodic limiter) has been later implemented in sev- 
eral large size and small size tokamaks. Particularly, an ergodic limiter has been 
installed in TEXT (Refs. [4-9]), in Tore Supra (Refs. [10-17]) and in a number 
of small-size tokamaks (in JFT-2M ( [18, 19]), CSTN-II ( [20,21]), HYBTOK-II 
( [22]), TBR-1 [23], TCABR ( [24]), and other small fusion devices ( [25,26]). 

The ergodic divertor implemented in the TEXTOR tokamak, the Dynamic Er- 
godic Divertor (DED) [27-29], has dynamical features in addition to the conven- 
tional ones. It permits the operation with a rotating magnetic perturbation field, 
which allows not only to broaden the heating footprints to avoid hot spots on the 
divertor plates but also to actively influence the plasma rotation and MHD activ- 
ities [30, 31]. At the same time the DED perturbation field has an influence on 
runaway electron dynamics [32,33], spontaneous density built-up [34], and reduc- 
tion of turbulence transport [35]. Particularly, the DED perturbation field has been 
used to control edge localized modes (ELM) in limiter H-mode plasmas [33], sim- 
ilarly to the recent experiments in the DIII-D tokamak to suppress ELM by the 
external resonant magnetic perturbations [37-39]. 

The energy and particle transport in the ergodic zone depends significantly on 
the structure of the perturbation magnetic field at the plasma edge. Numerous 
recent experiments on the TEXTOR-DED have demonstrated that the transport 
along the field lines in the presence of the resonant magnetic perturbations becomes 
predominant and prevails over the cross field neoclassical and anomalous diffusion 
[40]. The field lines even determine a fine structure of heat deposition patterns 
[41^3], which is mainly due to the reduction of the large scale turbulence level 
at the plasma edge during DED operation [35, 36]. Therefore the study of the 
magnetic field structure is an important first step in understanding the transport 
properties of the plasma edge modified by the resonant magnetic perturbations. 

Traditional methods to study the perturbation field structure and its statistical 
properties in ergodic divertors have been mainly based on numerical codes [44- 
47]. One of the approaches to study the formation of the ergodic zone of field 
lines is based on the direct field line tracing in a toroidal system (see, e.g., [46]). 
More in particular, to study the perturbation magnetic field structure in the ED 
of Tore Supra the MASTOC code has been developed. In its early version [44], 
the field lines equations are formulated in a Hamiltonian form using the intrinsic 
coordinates. The Fourier coefficients of the perturbation Hamiltonian are found 
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by direct integration along the unperturbed field lines using the Biot and S avail 
formula and the field line equations are integrated by the Runge-Kutta integration 
scheme. 

A new approach to study this problem with an application to the TEXTOR- 
DED has been developed in Refs. [48-51]. This approach consists of two steps. 
The first step is to formulate the equations of magnetic field lines in Hamiltonian 
form using the magnetic flux coordinates. For this purpose the analytical mod- 
els for the external coil configurations and the equilibrium magnetic field of the 
plasma are used [48,49,51]. In the second step the Hamiltonian equations are inte- 
grated using the computationally effective symplectic mapping methods developed 
in Refs. [52-54]). 

In this work we shall apply these analytical and mapping methods to study the 
perturbation magnetic field structure in the ED of Tore Supra. 

The paper consists of three sections. In Sect. 2 we consider a model for the 
ED perturbation coils, analytically calculate the perturbation magnetic field and 
study its poloidal and toroidal spectra, as well as its radial dependence. There we 
also study the poloidal spectrum of the perturbation magnetic field in magnetic 
flux coordinates, its relation with the corresponding spectrum in geometrical coor- 
dinates, and compare with the spectrum of the TEXTOR-DED perturbation field. 
The structure of the ergodic zone created by the perturbation field and the statistical 
properties of chaotic field lines are studied in Sect. 3 using different methods, rang- 
ing from the qualitative Chirikov criterium to the quantitative methods as Poincai^e 
sections, laminar plots and magnetic footprints. In addition, we calculate the radial 
profiles of the field line diffusion coefficients. The summary of obtained results is 
presented in the conclusive section. 

2 ED coil configuration, perturbation magnetic field and 
perturbation spectra 

Tore Supra is a large size tokamak with superconducting toroidal coils. Its major 
radius is Rq = 2.38 m and the minor radius a = 0.8 m. The ED coil configu- 
ration consists of six identical modules located on the low field side of the torus 
and equidistantly spaced along the toroidal angle ip [11,44]. Each module has a 
poloidal extension 27r/3 and a toroidal extention A(p ^ 2tt/1A. The mod- 

ules ai^e located at the minor radius Vc = 85 cm. 

We model the module with a coil winding shown in Fig. 1 where arrows indi- 
cate the cuiTcnt direction. The cuiTcnt flows from the feeder located at the begin- 
ning of the first section j = 1 of the inner side of the winding shown in Fig. la and 
returns through the outer side of the winding shown in Fig. lb. The minor radii 
of the inner and outer sides are rd = 84 cm = 86 cm, respectively ^ Below 
we calculate the magnetic field, created by the cuiTcnt flowing in this coil system 

'Furthermore, we use the quasitoroidal coordinate system {r,d,tp), which will be referred as 
geometrical coordinates. 
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(a) 



(b) 




Figure 1: Model scheme of one module of the ED coils: (a) the inner winding at 
Tci = 84 cm; (b) the outer winding at = 86 cm. 

by considering, in a first time, the case of inner and outer parts of the winding 
separately, and by summing those parts afterwards. 

Suppose that the modules are centered near the toroidal angles (pk = {k — 
l)Aip, Aip = 27r/6. The poloidal angle as a function of the toroidal angle of a 
point located in the j—th section in the k—th module (fc = 1, 2, . . . , 6) is given by 



{^) = ef'\^) + {j -1)56, j = l,2,...,N, 



(1) 



where A'' = 8 for the inner winding and = 6 for the outer winding, and 56 is the 
poloidal spacing shown in Fig. 1. The coordinates of the first coils on each module 
ai^e described by 



i^\^) = 6i{0) + a{^-ipk) 



for ipk 



^ < If < ipk + 



(2) 



where a is the slope of a coil creating a helical magnetic perturbation. 

The poloidal extension, A6, of a module shown in Fig. 1 can be expressed as 
a function of the poloidal position of the first coils, 0i(O), at the toroidal section 



A6 = 2\6i{0)\ + aA^p. 



(3) 



One should note that this model of coils is not fully equivalent to the Tore Supra 
coils. In the latter case the distance between sections of coils in each module are 
not equidistant along the poloidal angle 6. It slightly decreases with the distance 
from the equatorial plane ^ = [55]. 

2.1 Current density 

We describe the current, Ij, which flows in a coil section by 



if =/,cos(7^j) = (-ly+^ 



j = l,2,...,N, 



(4) 
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where is the current flowing in the coil, z = 1 for the inner part of the winding 
and i = 2 for its outer part. 

Below we shall consider only the long helical section coils since they create 
the magnetic field perturbations that are resonant with the magnetic field lines of 
the plasma. The vertical short sections of coils do not contribute to the resonant 
field, therefore they will not be taken into account. 

One can introduce the current density vector j (r, 9,ip) of the coil system as 

e/ N 6 TV 

fc=i j=i 

where e^*) = (eV, eg, eip) = (0, sin aoi, cos aoi) is a unit vector along the helical 
section of the coils, OQi = ard/Rc, Rc = Ro + fch = 1) 2). Here g^\'^) is 
a step function of the toroidal angle Lp which takes a non-zero value in the areas 
covered by coils, i.e.. 



' 1, for - ^ < 99 < + 



(6) 



0, elsewhere, 



Introducing the step function gi{6) solely depending on the poloidal angle 

g^{0) = 



"2 ^ " ^ 2 ' 



(7) 



0, elsewhere, 

the current density (5) can after some transformations be reduced to 

6 

rii) 



k=l 

00 

cos{moi2s-l)[{e-9o)-a{^-ipk)]}. (8) 



X 

s=—oo 



where 

One can show that the current density (8) can be expanded into a Fourier series, 

00 00 00 

Ur,9,ip) = 2e^'^ ^ ^ ^ ji^f) (r) cos [inO - + X^^i) , (9) 

m=— oon=— 00 s=l 

with the Fourier coefficients 
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Xs = mo(2s - 1)00, 



(10) 



where 

~M) ^ G^'^AifAe, ^ mold eA^AOj 

° (2vr)2 vrr,, (27r)2 ' ^''^ 

(,) _ sm{[m-mo{2s-l)]A6i/2) 
9mi [m-mo(2s-l)]A0i/2 ' ^ ^ 

^(s) _ siJ^ ([^ - "^o(2g - 1)q]Av?/2) 

[n-mo(2s- l)a]Av?/2 ' ^ ^ 

The toroidal mode number n takes values n = 6g, g = 0, ±1, ±2, .... As one 
can see from Eqs. (10) and (13), the biggest effect occurs when the ED coils are 
designed in such a way that the product rriQa is close to the toroidal mode no = 6, 
i.e., |?noa — no| ^ Aip/2. Then in the sum (9), the main contribution comes from 
the terms with the toroidal numbers n = (2s — l)no, (s = 1,2,...). Leaving in 
Eq. (9) only these terms we have 

oo oo 

ji{r,e,^) = 2(^'^ Y.^^rZl{r)xcos{me-{2s-l)noip + Xs), (14) 

m=— oo s=l 

where no = 6, i = 1, 2. 

2.2 Magnetic field in a cylindrical approximation 

First we consider the magnetic field created by the helical cuiTcnts (14). Using 
the procedure similar in Ref. [51] one can obtain the exact formula for the scalar 
potential <I>(r, 9, (p). Then the perturbation magnetic field is given by B(r, 6, ip) = 
V<l>(r, 6, ip). For r < r^, {i = 1, 2) one has 

$(r, e, ip) = $i(r, e, ip) + $2(r, 9, ip), 

oo 

Mr,e,p)= J2 Y.^^^nir)sm{me-np + xs), (15) 



m=— oo n 



where n = (2s — l)no, (s = 1, 2, . . . ) and 

$W (r) = fW (r) — 

mn\ I c i>m J mn\ ) ^ ' 

2fiomoldcos{aoi) 6A(pA9i 

^''^ 

where Irn{z) andi^m(z) are the modified Bessel functions {K'^{z) = dKm{z) / dz). 
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(a) (b) 




Figure 2: (a) Radial dependence of the function fmn{i^) for different poloidal 
modes m. (b) Relative deviation of the fmn{f) from the power law {r/rc)^: 
^fmn{r) = \ fmn{r) - (r/rj™ |//„„ (r). Parameters are = 85 cm, Rq = 238 
cm, n = 6. 



(i) J 

Here Be is the characteristic amplitude of external magnetic field strength . 
For the typical parameters of the ED of Tore Supra (rc = 0.85 m, = 22.5 kA) 
we have Be ~ 425 G. It is close to the value, B^ 405 G, of the con^esponding 
quantity for the TEXTOR-DED (at Id = 15 kA). 

The radial dependences of the perturbation field are described by functions 
fmn{r) which are shown in Fig. 2 for a several mode numbers m. For large 
mode number m {m > 4) the radial dependence is well described by the fol- 
lowing asymptotical formula fmn{r) ~ [r/rd)"^- Then the radial component of 
the magnetic field Br can be represented as 

Br{r, e, y,) = bP (r, 9, y,) + B^ (r, 9, 

oo oo / \ m— 1 



2.3 Toroidal corrections 

We consider the first order correction to the magnetic field perturbations due to 
toroidicity of the system. According to [56] this effect can be taken into ac- 
count by multiplying the scalar potential ^{r,9, ip) obtained in the cylindrical ap- 
proximation with the factor y'^^o7CR(7+^^"cos^ when the corrections of order 
{nrc/2Rc)"^ are small. In this case we have 



^' ^) = V a ^' ^) + 0, if)] , (18) 

\ Kq + r cos 9 

where the functions ^i{r,9,(p), {i = 1,2) are given by Eq. (15). We present 
the corresponding radial magnetic field as a sum of contributions corresponding to 



^We should note that the definition of Bc'^ as well as coefficients gin'' are slightly different from 
the corresponding ones given in Ref. [48, 51] 
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different toroidal modes, n, 



(19) 



Xs 



m=—oo 



where 



Bmnir,9) 



$(1) (r) + $(2) (r) 

mn\ I ' mn\ i 



d_ I Rq 
dr\jRo + r cos ( 

(\ m— 1 / \ m— 1 



rc2 



Rq 



r cos ( 



Rq + r cos \ 2m {Rq + r cos ( 



(20) 



The toroidal component of the vector potential of the perturbation field takes 
the form 

oo 

O'mni'^^ ^) cos ijnO — nip + Xmn 

m=—oo n 

amnir,0) = m~^rBmnir,0), 

(21) 

where e stands for the dimensionless perturbation parameter, defined as e = b'^'^ /Bq, 
Bq is the strength of the toroidal field. The dimensionless Fourier coefficients, 
o-mn{f, are then given by 



,ir,e) = agUr,6) + 6agUr,0), 



(22) 



where 



(0. 



r d 



Ro 



9.. 



mRo mdr \ ^ Rq + r cos \ 

(^) ( ^ \ .1 zz^ II 



' tuRq \ r, 



r cos 9 



Rq + r cos 9 \ 2m {Rq + r cos ( 



(23) 



Here (i = 1,2), 6 = bP/B^^^ = r.iA^a/ {rc2^9i)- The phase Xmn = Xs and 
toroidal mode number n = {2s — l)nQ. 

The angular- dependencies of the perturbation field Br{r, 6, ip) at the fixed val- 
ues of radial coordinate r and the toroidal angle ip = aie plotted in Fig. 3. 
The radial dependence of the perturbation field is determined by the poloidal mode 
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Figure 3: Radial component of the perturbation field Br', (a) contour plot in the 
(9, r) plane {^p = 0); (b) poloidal dependence at different radial positions r: 1 — 
r= 80 cm , 2 - r= 70 cm, 3 - r= 60 cm. 1^= 22.5 kA. 



spectra, gmi- According to the definition given by Eq. (12) the latter is localized 
near the central mode mc = (2s — l)mo = nniQ/riQ, (s = 1, 2, . . . ) and has a width 
Am = vr/A^j. Therefore, one expects that the perturbation field, Br"\r, 6, ip) of 
the n— th toroidal mode has the following radial dependence, 

4") (r, 9, if) oc (—) \ 7„ = ^ - 1. (24) 
\rcij no 

The power law of the radial decay of perturbation field has the lowest exponent, 
'yn=6 = niQ — 1, for the toroidal mode n = 6 . For the value 59 = 18°, one has 
mo = TT / 59 = IQ the exponent 7 = 9. For the next toroidal mode n = 18 we have 
7n=i8 = 3mo - 1 = 29. 

It should be pointed out that the power-law decay of the perturbation field (24) 
found here is different from the exponential decay law exp(— mo(rci — r)/rci) 
proposed for the ED in Tore Supra [44]. 

2.4 Determination of the spectra of magnetic perturbation Hmniip) 

Starting from the expressions (19) to (23) for the perturbation field Br{r, -d, ip) and 
the toroidal component of its vector potential, the field line structure in the Tore 
Supra ED edge can be studied following the approach described in [51]. In terms 
of magnetic flux coordinates {ippoi,tp = i/^tor , , ^) the magnetic field B can be 
presented in the Clebsch form [57-59], where {ippoW^) are poloidal and toridal 
fluxes, (i?, Lp) are poloidal and toroidal angles. In this coordinate system the field 
line equations have the Hamiltonian form 

dij) _ dippoi d-d _ dippoi ^25^ 
dip d-d dip dip ' 
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with the Hamiltonian 

^pol V') = e X] ■^"^'^(^) cos(?ni? -nLp + Xmn)- (26) 

m,n 

where q{Tp) is the safety factor of the equalibrium plasma, Hmn{'>P) the Fourier 
coefficients of the perturbation field conesponding to the poloidal mode numbers 
m and the toroidal mode n. The equations (25) are integrated using the computa- 
tionally efficient mapping method (see Refs. ( [52-54]). 

The equilibrium plasma of the Tore Supra tokamak is modelled by a plasma 
with nested, circular magnetic surfaces with an outward (Shafranov) shift due to 
effects of the plasma pressure and electric cun^ent, again completely similar to the 
equilibrium plasma model described in [51] (which is based on expressions for the 
Shafranov shift, the toroidal current density profile of a cylindrical plasma model 
and the representation of the safety factor due to toroidicity as a series of powers 
of the inverse aspect ratio given in [44], [60] and [48], respectively). The ripple of 
the magnetic field was not yet incorporated in the modelling presented here, thus 
ruling out the possibility of a direct comparison of the results with experimental 
evidence. 

In order to integrate the field line equations, the Fourier components Hmn{^) 
of the perturbation Hamiltonian have to be determined. According to Refs. [48,51, 
54] those are found by the following Fourier integral from the product of the major 
radius R = Rq + r cos 9 and the vector potential Aip of the perturbation field (21), 

27r 

'""^^^ JJ {2TryBoRl ^' 



where (r, 6) is the unperturbed field line on the given magnetic surface tjj and -& is 
the intrinsic poloidal angle which is different from the geometrical poloidal angle 
9. In the Eq. (27) it is supposed that (r, 6) is a function of the ip and ^: r = 
r{'il),'d),9 = 9{'il),'d). 

Using the Fourier series in Eq. (21) the integral can be reduced to 

m' 

/ f (30S \ ^^'^ 

hm'nir, 9) = am'nir, 9)il + ) • (28) 

As seen from Eq. (28) the toroidal spectra, n, of the perturbation field remain 
unchanged. However, the poloidal spectra, bmn{r,9), in the geometrical space. 
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(r, 9), changes drastically in the {ip, i}) space. First following Refs. [48,51,54] we 
qualitatively estimate the integral (28) in order to reveal the main features of the 
transformation of the spectra of magnetic perturbations. 

The integral (28) contains the fast oscillating phase = md — m'9{'d) and 
slowly varying amplitude bm/n{r,9). Therefore it can be estimated using asymp- 
totical methods [61]. As was shown in [51,54] for the perturbation field located on 
the low field side of the torus the Hmniip) is given by the following asymptotical 
formula 

E/ 2 \^/^ ( ^im' — m \ 
h^,n {r,e = 0) —— Ai --^ — — , (29) 

\\lz\m'J \ (|73|mV2)^/V 

where the coefficients 71 and 73 are the first and third derivatives of the geometrical 
angle 6 with respect to the intrinsic angle d taken on the low field side of the torus, 
= 0, respectively: 



de 



(30) 



We should note that 71 > 1, 73 < 0, and their particular values depend on the 
equilibrium plasma configuration. The sum (29) can be replaced by integration 
over m' and asymptotically estimated. This procedure gives 

HmnW ~ ll^Bbm'(^rn)n TT = 0) , (31) 

where m'{m) is determined by 

m-Xe(|73lW2)^^^ 
m [m) = , (32) 

71 



and B = y27r7|x^Ai {xc), and Xc ~ — 1 is the local maximum of the Airy func- 
tion. 

The obtained relation (31) describes the transformation law of the poloidal 
spectra of magnetic perturbations in a geometrical space to the ones in intrinsic 
coordinates. As was mentioned above the poloidal mode spectra a,„„ at the given 
toroidal mode n are localized near the central mode rric = nrriQ / no with a poloidal 
mode extension Am = -k/ /S.9i. Then from Eqs. (31), (32) it follows that the spec- 
trum, Hfnn, has the same form as Qmn, but its central mode m* is shifted to the 
higher number m* « mc7i and the width becomes larger Am* w Am7i, since 
7i > 1. These rigorously obtained results coincide with the corresponding ones 
obtained in [44, 55] by a qualitative analysis. 

However, the radial dependence of perturbation modes Hmn{tp) (31) has a 
power-law Hmni'4') oc r"^/"/^ in contrast to the exponential law exp {—■m{rc — r)/rc) 
supposed in [44,45] . 

The mode transformation is illustrated in Fig. 4 where the comparison of the 
asymptotical formula (31), (32) for Hmn{'4') (curve 2) with its numerically calcu- 
lated value (curve 3) from the integral (28) is shown. Pai^ameters B and Xc in the 
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Figure 4: Poloidal spectra of magnetic perturbation field at the magnetic surface 
q = 3: Curve 1 describes bm{r, 0), curve 2 corresponds to the asymptotical formula 
(31), (32) with the fitting parameters B = 1.4, Xc = —0.5, and curve 3 corresponds 
to the numerical calculations of the integral (28). The plasma parameters are the 
plasma minor radius a = 80 cm, the major radius Rq = 238 cm, the toroidal 
field Bt = 3.03 T, the plasma current /p=1.5 MA, the plasma Ppoi = 0.13. The 
ED current 1^= 22.5 kA, 66 = 19°, A(p = n/U. The quantities 71 = 1.85, 
73 = -3.124. 

asymptotical formula are considered as fitting parameters. Curve 1 corresponds to 
the spectra bm{r, 0) in the geometrical space. As seen from Fig. 4 the asymptotical 
formula describes well the transformation law of the poloidal spectra of magnetic 
perturbations. 

The asymptotical formula (31) resulting from the qualitative analysis of the 
transformation of perturbation modes in geometrical coordinates to the ones in 
flux coordinates is important, as it contributes to understanding the features of 
mode transformations in toroidal systems and clearly shows the differences be- 
tween the cases when the perturbation coils are located on the low-field-side (for 
the ED of Tore Supra) and on the high-field side (the TEXTOR-DED, for which the 
corresponding formula is given in Sect. 2.5), without having to resort to numerical 
computations. 

Fig. 5 shows the contour plot of the poloidal spectra of magnetic perturbations 
Hmn{ip) for the toroidal modes n = 6 and n = 18, respectively. The values 
Hmniip) the resonant surfaces (or ipmn, nq{'ipmn) = fn) lie at the white 
curve nq{p). One can see that the magnetic perturbation have a wide spectrum 
with the central modes m*{'>pmn) at the rational magnetic surface ijjmn, (qiipmn) = 
m/n) close to the corresponding resonant mode numbers m. It covers the rational 
magnetic surfaces located at the plasma edge between q > 2 and q = 3.5. The 
magnetic perturbation for the toroidal mode 77, = 18 is one order smaller than for 
the mode n = 6. 
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Figure 5: (a) Contour plot of poloidal spectrum of perturbation field, e| in 
the (m, p)— plane for the toroidal mode n = 6. White dashed curve corresponds 
resonant line {nq{p), p). (b) The same as in (a) but for the toroidal mode n = 18. 
Parameters are the same as in Fig. 4. 



The resonant components of Hrnn{'>Pmn) decay inwardly as Hmn{ipmn) cx 
i^mn , with nir ~ 7.9 for the toroidal mode n = 6. 

2.5 Comparison with the TEXTOR-DED 

Unlike the ED coils in Tore Supra, the perturbation coils of the DED of TEXTOR 
are located on the high-field-side (HFS) of the torus [48,49, 51, 62]. It consists of 
16 helical coils continuously winding around the torus with the poloidal extension, 
AO « 72°. By special arrangement of the current distribution over the coils one 
can create the perturbation field with the dominant toroidal modes, n = 4, n = 
2, and n = 1, which are called m : n=l2:4, 6:2, and 3:1 operational modes, 
respectively. The poloidal spectrum, \bmn{r,0)\, of the perturbation field in the 
geometrical space has a wide spectrum with the central mode rric = nruQ /4 5n 
and the width Am ^ ir/ A9, where niQ f» 20. It is shown as curve 1 in Fig. 6 for 
the toroidal mode n = 4. The radial decay of the modes \brnn\ has a power law 
\bmn{'f,0)\ oc (r/rc)™, and con^espondingly the full radial perturbation magnetic 
field, Br oc {r /r^)'^^ , with the exponent, 7„ mo — 1 = 19, which is twice larger 
than the corresponding exponent 7^ for the ED of Tore Supra (see Eq. (24)). It 
means that the perturbation field in the TEXTOR-DED does not penetrate into the 
plasma much deeper than in the case of the ED of Tore Supra. 

On the other hand the transformation from geometrical to intrinsic coordinates 
for the TEXTOR-DED case, modifies the perturbation field spectrum, €\Hran\, in a 
completely different way. An example of this spectrum is shown in Fig. 6 as curve 
2. The central mode number niQ of the poloidal spectrum in the geometrical space 
is about mo ~ 20 and the width Am k, 10. In intrinsic coordinate the central 
mode is shifted to ml k, 12, and the width Am* « 6. 
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Figure 6: Typical poloidal spectrum of perturbation field e\Hrnn\ at the mag- 
netic surfaces g = 3 in the TEXTOR-DED: curve 1 describes bm{r,TT), curve 
2 corresponds to numerically calculated Hmn{i^)- Parameters of the plasma are 
Bt = 1.9 T, Ip = 382 kA, (i^oi = 0.5, a = 44.7 cm, Rq = 175 cm, the DED 
current, 1(^=15 kA. 

As was shown in Refs. [48,49,51] this spectrum has the following asymptotics 

m (m) = , (33) 

Pi 

where (3i and /^s are the first and the third derivatives of 9 with respect to d taken 
on the high field side ?9 = vr [51,54]. We note that /3i < 1. In this case the central 
mode ?Ti* is shifted to the lower mode number m* = rricPi, and the width of the 
spectra is decreased. Am* w Am(3i . 

Thus the analytical analysis made in Sees. 2.4 and 2.5 show that the effect 
of the perturbation significantly changes when the coils are moved from one side 
of the torus to the opposite one. If the TEXTOR-DED coils were located on the 
low-field-side with the same coil design the effect of the perturbation would be 
drastically increased: the group of modes of bm centered neai^ rric ~ 20 would be 
transfered to the mode group with the central mode number m* w ^rric ~ 36 and 
the width Am* w 7Am ^ 18 which however are not resonant at the plasma edge. 
The resonant modes of Hmn with m w 10 — 14 would be formed from the group 
of modes of 6^ with m' = 5 — 8. The radial decay of these modes is much weaker, 
i.e., Hmn oc r''', (7 = 5 — 8), and they penetrate much deeper into the plasma. 
In the case of Tore Supra the move of the ED coils to the opposite side would 
significantly weaken the effect of perturbation. The resonant modes Hmn with 
m = 14 — 20 will be formed from the group of modes of b'^ with m' w 25 — 33, 
which has a much stronger radial decay Hmn oc r^, (7 = 25 — 33). 

In reality the coils for the Tore Supra ED as well as the TEXTOR-DED are 
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designed to have a group of modes near the central mode perturbation, m*, which 
are resonant to field lines near the q = 3 surface of the plasma edge. 

In addition, from Eqs. (31) and (33) the following roughly determined radial 
decay laws for the perturbation modes can be derived: H-mnW oc r™/^i, (71 > 1) 
for the ED of Tore Supra and H^^iip) oc r'"/^! , < 1) for the TEXTOR-DED, 
i.e., in the first case at a given poloidal mode in, the radially inward decay of 
the perturbation is much weaker than in the second case. This is most probably 
related to the location of perturbation coils on the low-field-side and the high-field- 
side, respectively, although differences in design between the perturbation coils in 
Tore Supra ED and TEXTOR-DED can not completely be ruled out as a possible 
explanation. 

3 Structure and statistical properties of the ergodic zone 

3.1 Qualitative estimates of the formation of the ergodic zone at the 
plasma edge 

Similai^ly to the analysis which has been made for TEXTOR-DED [51], a qualita- 
tive picture of the formation of the ergodic zone in Tore Supra ED can be obtained 
by considering the Chirikov parameter [63,64], which is determined by the pertur- 
bation spectrum and the q(/o)-profile. 

This Chirikov parameter, achir, as a function of the mean radius of the neigh- 
boring resonant magnetic surfaces, p = {pm+i,n + P-mn)/'^, is shown in Fig. 7 for 
the three different perturbation currents: Id = 4.5 kA (curve 1), Id = 9 (curve 2), 
and Id = 22.5 kA (curve 3). As seen from Fig. 7 the Chirikov parameter, acur^ 
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Figure 7: Chirikov pai^ameter achir vs the radial coordinate p for the different 
perturbation currents: curve 1 corresponds Id = 4.5 kA, curve 2 — Id = 9 kA, and 
curve 3 — Id = 22.5 kA. The plasma pai^ameters ai^e the same as in Fig. 4. 

grows linearly with the radius, p, and exceeds the unity (which is the criterium 
for the overlapping of magnetic islands) at p > pi 66 cm for the perturbation 
current Id = 4.5 kA, p > /32 ~ 64 cm for Id = 9 kA, and p > ps « 61 cm 
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for Ifi = 22.5 kA. Therefore, the field lines ai^e chaotic in the regions, p > pi, 
{i = 1, 2, 3), at the corresponding perturbation currents. Id = 4.5 kA, Id = 9 kA, 
and Id = 22.5 kA, respectively. 

3.2 Structure of the ergodic zone 

A more precise picture of the onset of chaotic field lines and the formation of the 
ergodic zone can be obtained by Poincare sections and laminar plots of field lines. 
Below we construct these plots employing the earlier mentioned mapping method 
of integration of field line equations. Furthermore, the integration step, A(p, of the 
mapping along the toroidal angle ip, is chosen equal to In/ns, with n = 6 and 
s = 4. 

In Figs. 8a-c Poincare sections of field lines in the poloidal plane 99 = are 
plotted for the three different perturbation currents. Id', (a) Id = 4.5 kA, (b) Id = 9 
kA, and (c) Id = 22.5 kA for the plasma and ED parameters corresponding to the 
spectrum of magnetic perturbations shown in Fig. 4. We plotted Poincare sections 
in the (p, -d) plane which clearly reveal the positions of the magnetic surfaces. (In 
this plane the unperturbed field lines are presented by straight lines, p =const). 

We can see from Fig. 8a that even at the lower level perturbation current 1^=4.5 
kA one obtains a rather highly-developed ergodic zone in the radial region p > 
Pc « 63 cm with open field lines to the plasma wall. It is formed by the overlap 
of nine primary neighboring magnetic islands with the poloidal mode numbers 
from m = 11 to m = 19 and the toroidal mode n = 6. There are only two 
remnants of KAM-stability islands near the boundary of the ergodic zone with the 
last intact magnetic surface pc ~ 63 cm. In the radial direction the ergodization 
level increases and the remnants of KAM islands drastically shrink in size. The 
region of field lines with short wall to wall connection lengths (white areas in the 
Poincare plots) grows when one approaches the plasma edge p = 80 cm. This area 
is known as a laminar zone. 

With increasing perturbation current Id the ergodic zone grows radially inward. 
At the same time the radial width of the laminar zone increases (see Fig. 8b,c). At 
the maximal perturbation current 1^=22.5 kA the magnetic island m = 10, n = 6 
is overlapped to the ergodic zone and the radial extension of the laminar zone 
becomes Api 12 cm. (A more precise way of estimating the laminar zone 
boundary will be given in Sect. 3.5.1 using the radial dependence of the field line 
diffusion coefficients). 

The magnetic field lines in the ergodic divertor are a typical example of an open 
chaotic scattering system and it is characterized by a certain (fractal) structure, 
especially, in the laminar zone (see, e.g., [51,65]). The structure of the latter can 
be more clearly visualized by the so-called laminar plots, i.e., the contour plots 
of field line wall to wall connection lengths. Such a plot of connection lengths 
measured in poloidal turns is shown in Figs. 9a,b in the {-d, p) plane at the maximal 
perturbation current 7^=22.5 kA: (a) shows the plot in the whole poloidal section 
and (b) shows the expanded view of the rectangular ai^ea at the low-field-side (LFS) 
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Figure 8: Poincare sections of field lines in the p) plane at tiie tiiree different 
perturbation cun^ents: (a) — Id = 4.5 kA, (b) — = 9 kA, and (c) — Id = 22.5 
kA. The toroidal modes n = 18 beside the main mode n = 6 is taken into account 
into account. The plasma parameters are the same as in Fig. 4. 



shown in (a). The corresponding Poincare section was shown in Fig. 8c. 

As seen from this figure the laminar zone forms a regular lattice-like structure 
in the large poloidal extended area, except the region on the LFS where field lines 
are distorted because of the magnetic perturbation which is localized on this side 
of the torus. One notices that the light blue "cells" in the laminar plot which cor- 
respond to the regions of field lines with wall to wall connection length of one 
poloidal turn are regularly located along the poloidal and radial directions. They 
are positioned between the resonant radii, Pmn- Their number along the poloidal 
direction at the given radius, p, coincides with the poloidal mode number, m, of 
the resonant radius, Pmn, located below the "cells". The sizes of the "cells" grow 
with increasing radius p. 

The expanded view of the laminar plot at the LFS shows also the so-called 
private flux zone (dark blue areas) of field lines which do not enter into plasma. 
Field lines coming from the inner regions of the plasma are connected to the wall 
along the elongated stripes (they aix seen as black stripes). 
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Figure 9: Laminar plots of field lines in the (iD, p) plane at the maximal perturbation 
current Id = 22.5 kA corresponding to Fig. 8 c: (a) on the whole plane; (b) shows 
the expanded view of the rectangular area in (a). 



3.3 Comparison with the TEXTOR-DED 

As seen from the spectra of the magnetic perturbation, Hmn shown in Fig. 6 the 
ergodic zone at the plasma edge in the TEXTOR-DED is formed by the interaction 
of only a few poloidal modes, 11 < ?n < 14, (12:4 operational mode) in contrast 
to the ED of Tore Supra (11 < m < 20). On the other hand, in TEXTOR-DED the 
magnetic perturbation radial decay has a power-law dependence with the exponent 
7 twice as large as in the ED of Tore Supra. Because of that one obtains in normal 
TEXTOR discharges an ergodized zone of field lines which are weakly connected 
to the wall. In order to increase the magnetic flux due to the perturbation field the 
plasma column is shifted to the HFS from the center of the torus. 

An example of the ergodic zone when the plasma center is shifted over 3 cm 
from the center of the torus into the direction of the perturbation coils is shown in 
Fig. 10 where the Poincare section and the laminar plot of field lines are displayed. 
The plasma parameters correspond to the discharge #93100: Ip =382 kA, Bt = 
1.9 T, Ppoi = 0.5. 

As seen from Fig. 10 in this shot the region of the stochastic field lines mostly 
consists of laminar- zone with short wall to wall connection lengths. The width of 
the laminar zone is radially increased in expense of the ergodic zone of field lines 
with long connection lengths. 

3.4 Magnetic footprints 

Heat and particle deposition patterns on the divertor target plates are closely linked 
to the magnetic footprints on the corresponding plates. Since the geometry of 
divertor plates of Tore Supra is rather complicated, and in order to simplify the 
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Figure 10: Poincai^e section (a) and laminar plots (b) of field lines in the (??, p) 
plane for the TEXTOR-DED. The plasma parameters are the same as in Fig. 6. 
The perturbation current = 7 kA. 



problem we study in a first approach the structure of magnetic footprints on the 
imaginary divertor plate at the minor radius = 80 cm. We present them as 
contour plots of the wall to wall connection lengths measured in poloidal turns 
on the surface of rd = 80 cm. The procedure of corresponding calculations is 
described in Refs. [51]. 

Such a plot of magnetic footprints is given in Fig. 11 for the same plasma and 
ED current parameters as in Fig. 8 c. Fig. 11a describes the magnetic footprints 
on the whole toroidal surface = 80 cm, and (b) shows an expanded view of the 
region < ip < 60°, —120° < 6 < 120°. As seen from these figures the footprints 
consists of n = 6 helical stripes con^esponding to the toroidal mode number n as 
in the case of the TEXTOR-DED. Field lines coming from the inner plasma region 
hit the divertor plate along these helical stripes. Dark blue areas correspond to the 
private flux zone from which field lines cannot penetrate into the plasma. 

If, in a more detailed approach, we take the exact geometry of the Tore Supra 
divertor plates into account, we need to adjust for two different aspects. Firstly, 
we need to modify the criterium used to determine if field lines hit the wall after 
each map step. Secondly, the radius r of the field line mapping starting positions 
on the ED module surfaces, will be a function of the toroidal and poloidal angle. A 
footprint plot in the {4>,0) - plane taking those calculation procedure modifications 
into account, has been made for ED module 1 as a whole (shown in Fig. 11c) 
and, with higher resolution, for neutralizer 4 of the same module (Fig. lid). On 
Fig. 1 Ic the edges of the neutralizer plates have been indicated and the neutralizer 
plates have been numbered. On Fig. 1 Id we notice the pattern of helical stripes of 
higher connection length. The bended form of those stripes is due to the projection 
on the (j),6 - plane of the spatial structure of the neutralizer surface. It is clearly 
visible that, analogously to what was found for TEXTOR-DED [65], the helical 
stripes exhibit a fractal structure, i.e. self-similarity at different spatial scales. 
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Figure 1 1 : Contour plots of wall to wall connection lengths (magnetic footprints) of 
field lines in the imaginary divertor plate = 80 cm The perturbation cun^ent = 
22.5 corresponding to Fig. 8 c: (a) on the whole plane; (b) shows the expanded view 
of the region < ip < 60°, -120° < 9 < 120°. (c) footprint plots on the surface 
of real divertor plate for ED module 1 ; (d) an expanded view of the above picture 
on the neutralizer 4. 
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3.5 Statistical properties of field lines 

In order to estimate the level of radial heat and particles transport induced by the 
chaotic field lines one should estimate the statistical characteristics of those field 
lines. More precisely, one should study the radial diffusion coefficients of those 
field lines and the Kolmogorov lengths which characterize the degree of divergence 
of neighboring field lines. Below we calculate these statistical chai^acteristics and 
compare them with their quasilinear estimations. For this we will follow the corre- 
sponding procedure described in Ref. [51]. 

3.5.1 Diffusion of field lines 

Consider the second moment of the radial displacement of field lines, 

1 ^ 

1=1 

where / = Ro^f is a length of field lines. Averaging in (34) is performed over a 
number N of field lines uniformly distributed at the initial ip = plane on the 
given magnetic of radius p. The local radial diffusion coefficient, DFL,is defined 
as 

DfUp) = (t\1)/21, 
at the initial linear growth regime ct^(/), as illustrated in Fig. 12. The diffusion 
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Figure 12: Determination of a local diffusion coefficient, DpL, from the linear part 
of the dependence of the second moment of the radial displacement, cT^(p), of field 
lines on their length: curve 1 corresponds to (T^(p)/2/s, vs //Z^; curve 2 describes 
the lineal" curve DflI/Is- Here Is is a length of a field line per one map step. 

coefficient determined in this way will be a local function of radius p. 

Diffusion of field lines in the highly developed ergodized zone can also be 
estimated in the frame of quasilinear theory. One can show (see, [66, 67]) that for 
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the Hamiltonian system of field lines mentioned in 2.4 and extensively described 
in [51] the distribution function, f{ijj), of field lines along the toroidal flux, Tp is 
described by the Fokker-Planck equation 



where the diffusion coefficient Dm{iP) is given by the quasilinear formula 



The diffusion coefficient, Dq (p) , along the radial coordinate, p, is related to Dm [ifj] 
as 



The radial dependencies of the diffusion coefficients Dfl{p), numerically 
calculated from the field line equations and the corresponding quasilinear diffusion 
coefficients, Dq{p) for the three different perturbation currents. Id, are displayed 
in Fig. 13: 1^ = 4.5 kA (curves 1 and 4), 1^ = 9.0 kA (curves 2 and 5), and 
Id = 22.5 kA (curves 3 and 6), (solid curves 1—3 describe DpL{p) and dashed 
curves 4—6 — Dq{p)). The corresponding radial profiles of the averaged field line 
connection lengths are shown in Fig. 14. Again, the averaging was performed over 
a number N of field lines uniformly distributed at the initial = plane on the 
given magnetic of radius p. 

One sees from Fig. 13 that the numerically calculated diffusion coefficients, 
DpL, steadily grow in certain radial regions, pmin < P < Pi- Close examination 
of the corresponding Poincare sections given in Figs. 8 and the averaged connec- 
tion lenghts shown in Fig. 14 indicates that these radial regions correspond to the 
highly developed ergodic zone of field lines possessing the long connection lenghts 
with the its lower boundary at pmin and its boundary with the laminar- zone at pi. 
However, starting from pi, the numerical diffusion coefficients. Dpi, stop to grow 
and abruptly go down, while the quasilinear diffusion coefficients, Dq{p), still 
continue to increase. The radius, pi, can be called as the lower boundary of the 
laminar zone. As seen from Fig. 14 in the laminar zone p > pi the field lines have 
rather short connection lenghts. 

It is remarkable that for the smaller perturbation currents. Id = 4.5, 9 kA, 
the numerical diffusion coefficients, Dpi, perfectly follow the quasilinear formula 
(37) in the ergodic zones of field lines with long connection lengths. However, they 
sharply deviate in the laminar zone with the short connection lengths. This con- 
trasts with the radial profiles of diffusion coefficients presented in [51] for the case 

^Note that the factor 1/2 in the definition of the quasilinear diffusion coefficient given in Ref. [66] 
is absent. This factor is included into the diffusion equation. 

''The formulas (36), (37) for the diffusion coefficient are given in a singular form. Their ac- 
tual values should be taken at the resonant values of p = pmn (or ^ = ^mn)'- DAii^mn) = 




(35) 




(36) 



Dq{p) = {Rl/p^)DM{i^) 



(37) 



{l/2)TVe^q{lpmn)m'^\Hmn{lpmn)\'^ ■ 



22 



X 10"' 




P [cm] 



Figure 13: Radial dependencies of local field line diffusion coefficients, DpL, of 
field lines at the thiee different perturbation cuiTcnts: curves 1 and 4 con^espond 
to Id = 4.5 kA, curves 2 and 5-/^ = 9 kA, and curves 3 and 6 — = 22.5 
kA. Solid curves 1—3 describe the numerical calculations and dashed curves 4—6 
describe the quasilinear values, Dq. 

of the TEXTOR-DED 12:4 operational mode, where the correspondence between 
numerical and quasilinear diffusion coefficients in the ergodic zone is significantly 
worse. This difference might be explained by noticing that for the latter case, the 
applicability of the quasilinear theory is seriously hampered because the ergodic 
zone is formed by the overlap of only a few neighboring magnetic islands. In the 
ED Tore Supra case, on the contrary, a higher number of resonances, located closer 
to each other and thus showing a larger degree of overlap, result in a much better 
developed ergodic zone. 

On the other hand, for the maximum perturbation current 1^ = 22.5 kA the sit- 
uation is different. The difference between the numerical DpL and the quasilinear 
diffusion coefficient monotonically increases with the radius already in the ergodic 
zone, and the numerical DpL{p) exceeds Dq{p). The reason for such behavior 
of the Dpl{p) is not clear yet and has to be investigated. A first indication might 
be given by the Poincare section of Fig. 8 (c), where the ergodic character of the 
radial region just inward of pi seems to be less pronounced than in the cases with 
the smaller perturbation currents. This suggests the existence of some additional 
averaged outward drift, due to a deeper reaching and dominating laminar zone. 

In the case of the TEXTOR-DED the comparison of the numerical Dfl{p) 
with the quasilinear diffusion coefficients, Dq{p), is not quite justified because 
the ergodic zone in the DED is formed by the interaction of only three or four 
magnetic islands. The quasilinear theory of diffusion can hardly be applied to such 
an ergodic zone. Particularly, for the case shown in Fig. 10 the numerical diffusion 
coefficient Dfl{p) reaches the maximal value 3.7 x 10~^ m^/m at p = 42.95 
cm, while its quasilinear value is 4.74 x 10~^ m^/m at the same radius, i.e., the 
quasilinear formula overestimates the diffusion coefficient by a factor 10. 
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Figure 14: Radial dependence of the averaged field line connection length, mea- 
sured in poloidal turns, at the three different perturbation currents: curve 1 corre- 
sponds to Id = 4.5 kA, curve 2 — Id = 9 kA, and curve 3 — Id = 22.5 kA. The 
region nearest to the plasma edge has been enlarged in the upper right corner. 



3.5.2 Kolmogorov lengths 

The Kolmogorov time is a statistical characteristic of dynamical chaotic systems 
which characterizes the time of loss of information on initial state of the system. 
For chaotic field lines it corresponds to the characteristic decay length of the corre- 
lation of neighboring field lines, and called the Kolmogorov length. The calculation 
of the latter is based on the estimations of the Lyapunov exponents. 

The numerical procedure of finding Kolmogorov lengths, Lk, for the chaotic 
field lines in the ergodic divertor is described in [51,68]. It is based on the determi- 
nation of the local Lyapunov exponents, <j{p), describing the degree of exponential 
divergency of chaotic field lines launched from a given magnetic surface of radius 
p. The Kolmogorov length, Lk, is inverse proportional to a{p) averaged over the 
given magnetic surface, i.e., Lk = l/a{p). 

The analytical formula for the Kolmogorov length, Lk, has been proposed by 
Ghendrih et al. [?], which expresses Lkq, through the Chirikov parameter, 

Lkq = nqRo {ixachv,-l'i)~^"' • (38) 

The numerically calculated radial profiles of Lk for the chaotic field lines in 
the ED of Tore Supra ai^e shown in Fig. 15a for the three perturbation currents: 
curve 1 corresponds to Id = 4.5 kA, curve 2 — Id = 9 kA and curve 3 — 22.5 kA. 
The corresponding profiles of the analytical formula (38) are shown in Fig. 15b. 
The numerically calculated Kolmogorov length, Lk, is of the order of the con- 
nection length, Lc = 2TTqRQ, which equal 44.9 m at the q = 3 magnetic surface. 
However, a comparison of Fig. 15a and b shows that the numerical Kolmogorov 
length exceeds the the analytical estimation (38) by factors between five or ten. 
This disperancy is probably related to the approximate character of the formula 
(38). However, the numerical Kolmogorov length is more close to the formula 

—4/3 

Lkq ~ 2-nqRoa(j^^^. given in [45]. 
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Figure 15: Radial profiles of Kolmogorov lengths, Lk, of field lines: (a) calculated 
numerically; (b) according to the quasilinear formula. Curves 1 to = 4.5 kA, 
curves 2 and 5-/^ = 9 kA, and curves 3 and 6 — 1^ = 22.5 kA. 

4 Conclusion 

In summary we have studied the magnetic field structure of the ED of Tore Supra 
through the application of the asymptotical and mapping methods developed for 
the similar- study of the DED of TEXTOR. To this end, an analytical model for the 
ergodic divertor of Tore Supra that would recover the main features of its magnetic 
field structure has been proposed. The fine details of the magnetic field pertur- 
bations, its poloidal and toroidal spectra in the geometrical coordinate system as 
well as in straight-field-line coordinates, the structure of the ergodic and laminar 
zone, and the field line diffusion coefficients were investigated. The use of standard 
code methods similar to the MASTOC code to study these detailed aspects of the 
magnetic field in the ED would be significantly complicated by the requirement of 
much higher computer power or computation times. 

The decay of the magnetic field perturbations for the ED in Tore Supra was 
found to exhibit a power-law dependency on the radius, in contrast to an originally 
proposed exponential decay law. On the other hand, rigourously obtained expres- 
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sions for the width and the central mode of the perturbation magnetic field poloidal 
spectrum do coincide with the corresponding qualitatively obtained original re- 
sults. An asymptotical formula describing the transformation law of the poloidal 
spectra was also obtained. Compai^ed to TEXTOR-DED, where the ergodic zone 
at the plasma edge is formed by the interaction of only a few poloidal modes, the 
ED has a wide spectrum and for a given poloidal mode the radially inward decay 
of the perturbation was found to be much weaker. 

For the Tore Supra ED, Poincai^e sections, laminai" plots of chaotic field lines, 
as well as the calculation of field line diffusion coefficients showed a highly- 
developed ergodic zone at the plasma edge even for relatively low perturbation 
currents. Similar to the case of the TEXTOR-DED, the obtained footprint plots ex- 
hibited helical stripes coiTcsponding to the toroidal mode number and with a fractal 
structure. It was found that in the ergodic zone with long wall to wall connection 
lengths the numerical diffusion coefficients of field lines are well described by the 
quasilinear formula, while in the laminar zone with short wall to wall connection 
lengths they sharply deviate. However, for the maximum perturbation cuiTcnt, a 
different situation occured. Already in the ergodic zone, the difference between 
the numerical and the quasilinear diffusion coefficient, the former exceeding the 
lather, monotonically increased with the radius. We have also numerically calcu- 
lated the Kolmogorov lengths of field lines. They have an order of the connection 
lenghts. However, the numerically determined Kolmogorov lenghts exceeds the 
quasilinear- prediction by factors between five and ten. 
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